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Abstract 

We describe mixing scalar particles and Majorana fermions using Closed-Time- 
Path methods. From the Kadanoff-Baym equations, we obtain the charge asym- 
metry, that is generated from decays and inverse decays of the mixing particles. 
Within one single formalism, we thereby treat Leptogenesis from oscillations and 
recover as well the standard results for the asymmetry in Resonant Leptogene- 
sis, which apply when the oscillation frequency is much larger than the decay rate. 
Analytic solutions for two mixing neutral particles in a constant-temperature back- 
ground illustrate our results qualitatively. We also perform the modification of the 
kinetic equations that is necessary in order to take account of the expansion of the 
Universe and the washout of the asymmetry. 
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1 Introduction 



There are several different methods of calculating the lepton asymmetry, that may be 
generated in the Early Universe. This variety of approaches is a feature that Leptogen- 
esis [H has in common with many other topics in Particle Theory. 

The approach that is perhaps most commonly applied is to calculate the decay asym- 
metry of individual right handed, singlet Majorana neutrinos from S-matrix ele- 
ments These asymmetries, obtained by Quantum Field Theoretical methods, are 
subsequently substituted into classical Boltzmann equations, that describe the macro- 
scopic non-equilibrium dynamics and the evolution of the asymmetry. By LSZ-reduction, 
S'-matrix elements can be obtained from time ordered n-point Green functions. We then 
enjoy the usual simplifications of time-ordered perturbation theory, i.e. we can calculate 
the n-point functions using the standard diagrammatic Feynman rules. The price to pay 
for this simplification is cutting the direct link between the non-equilibrium dynamics 
and the generation of the CP-asymmetry at the quantum level. A particular problem 
turns out to be due to the instability of the singlet neutrinos, which causes that the 
S'-matrix, that involves the CP- asymmetric decays and inverse decays of these singlets, 
is not unitary, or rather, that it does not correspond to an S'-matrix in the proper sense. 
As a consequence, the unitary evolution of the non-equilibrium system must be care- 
fully re-implemented, which is performed in practise through the non-trivial method of 
real intermediate state subtraction jsj . Without real intermediate state subtraction, one 
would predict from substituting the rates from the S'-matrix elements into the network 
of Boltzmann equations the presence of a lepton asymmetry even in thermal equilibrium, 
which would be an unacceptable violation of the CPT theorem. 

For Resonant Leptogenesis it is assumed that the mass difference AM 

of two of the singlet neutrinos is small compared to their average mass M, or, more 
precisely, the energy difference is small compared to the average energy. Decays and 
inverse decays, that induce off-diagonal correlations of the singlet neutrinos in their 
mass eigenbasis, then lead to a relatively large mixing, which enhances CP- violation. In 
the S'-matrix approach, the mixing appears as a wave function correction to the singlet 
neutrino propagator in the decay diagrams BliMi. The mixing that is calculated in 
this way is time-independent, which excludes the possibility of incorporating the effect 
of oscillations. 

Alternatively, one can avoid the subtleties concerning unitarity by computing the 
real-time evolution of the quantum state from the outset. This can either be achieved by 
using the canonical Hamiltonian approach or functional. Closed- Time-Path (CTP) |l4j - 



ITj methods. This may appear as a complicated task, but one should keep in mind that 



the CP-asymmetry arises in the S'-matrix language at the one-loop level. In a real time 
approach, it turns out that in a diagrammatic representation, we have to go to two-loop 
level. However, the leading CP-violating effects can be extracted by imposing some of 
the internal particles to be on-shell, which may be considered as a generalisation of the 
Cutkosky rules 18j. In the end, the integrals that are relevant for the leading order 



prediction of the CP asymmetry are identical to those encountered in the S'-matrix 
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approach, but with the advantage that an overcounting of particular processes in the 
Boltzmann equations, that would violate unitarity, is avoided by construction. 

The canonical approach makes use of the property of the Hamiltonian as the time 
evolution operator. A blueprint of this canonical approach, is the theory of neutrino os- 
cillations in Refs. 19l-l23|. It has been noticed that when supplemented by CP- violation, 
the oscillations may lead to the creation of a lepton asymmetry 2J-|26| . The possibility 



of describing mixing oscillations is one key advantage of the canonical approach when 
compared to using S'-matrix elements. On the other hand, without a set of diagrammatic 
rules, the perturbation expansion is rather cumbersome. This is why the remaining im- 
portant CP-violating contribution, that is sometimes referred to as direct CP-violation 
{e.g. in Ref. [gI), has not been calculated in the Hamiltonian approach. In the 5- matrix 
approach, this contribution arises from the penguin-shaped vertex diagram. Besides, the 
decays of the singlet neutrinos into Standard Model leptons are described in practise 
by appealing again to the S'-matrix formalism. Therefore, one may state that within 
the Hamiltonian approach, no entirely closed picture for the emergence of the lepton 
asymmetry has yet been developed. 

In principle, it should be possible to derive a set of diagrammatic rules that simplifies 
the calculations within the Hamiltonian formalism. However, such a set of rules is 
readily available. It arises in an elegant and rather intuitive manner from the functional 



formulation of the OTP formalism |l4Nl7l|. In the OTP approach, the Hamiltonian 



time evolution of the density matrix is replaced by Schwinger-Dyson equations for the 
two-point Green functions. These equations have four matrix components from the 
two indices that denote the two branches of the CTP. The Kadanoff-Baym equations are 
those of the Schwinger-Dyson equations, that encompass the macroscopic evolution of the 
system. Recently, this approach has proved successful in order to describe Leptogenesis 
within a single theoretical framework 27H37l|. Both, CP- violation from wave- function 
and vertex corrections are incorporated. Unitarity issues are resolved and an accurate 
account of all quantum-statistical effects on the asymmetry is made. Moreover, the 
formulation in terms of Green functions bears the potential of incorporating corrections 
from Thermal Field Theory within this formalism. For the CP-conserving production 
rate of singlet neutrinos, several phenomenologically relevant thermal corrections have 
recently been calculated in Refs. 38N44j . When it comes to the rate of CP- violation, 
an aspect that is not dealt with in Refs. 27H37| is the possibility of flavour oscillations 
between the singlet neutrinos. The present work addresses this shortcoming. 
While in Ref. 
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the wave function correction is calculated using the perturbative 
loop expansion, the present work captures the same effect by solving the Kadanoff-Baym 
equations for the singlet neutrino propagator directly. In the limit where the mass dif- 
ference is much larger than the decay rate Fd of the singlets, AM ^ Fd, the standard 
perturbative result is reproduced. When AM ~ Fd or smaller, it becomes important 
that the solution to the Kadanoff-Baym equations corresponds to a resummation of all 
wave-function insertions. This is a generalisation of the familiar procedure for time- 
independent situations, which can be performed by the summation of a geometric series. 
The result obtained in the present work from the Kadanoff-Baym equations encompasses 
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both, time-independent and oscillating contributions to the mixing of the singlet neutri- 
nos. We note that for a model of the asymmetric decay of a scalar field, it was clarified 



in Refs. [45|, |46| that the Hamiltonian approach incorporates indeed the wave-function 
contribution to the CP-asymmetry from decays in the vacuum as it is usually calcu- 
lated from S'-matrix elements. The present work goes beyond this proof of principle 
not only because it also treats the decays of Majorana fermions, but also because it for- 
mulates the evolution of the asymmetry in an interacting, finite-density system within 
the single framework of the CTP formalism. Our results may therefore be employed for 
phenomenological studies of Resonant Leptogenesis in the Early Universe. 

In order to further describe the context of this present work, we note that for compu- 
tations of the asymmetry from S'-matrix elements, it has been realised that the resonant 
enhancement factor 1/AM should be regulated by the finite width of the singlet neu- 
trinos 0, E lO, 47 1 . For solutions to the Kadanoff-Baym equations in Wigner spacejl] it 



turns out that the non-equilibrium distribution functions are not multiplied by a finite- 
width, Breit- Wigner distribution, but rather by a zero width, Dirac (5-function. The 
observation of this somewhat curious property has lead to the speculation that these 
zero-width solutions arise due to an incorrect treatment of pinch singularities that must 



be addressed by a non-perturbative resummation [48|, |49|, that is yet to be specified. 
This may be perceived as an obstacle for the formulation of an effective theory of Res- 
onant Leptogenesis within the CTP framework. However, it was recently demonstrated 
analytically and numerically, that the vanishing width is a correct property of the out- 
of-equilibrium Wigner functions. Within loop diagrams, the zero width is corrected for 
by a resummation of derivative operators, which effectively shifts the argument of the 6- 



function from the real axis into the imaginary direction of the complex plane [50| . In the 
Kadanoff-Baym equations, this effect may be of relevance when the tree-level scattering 
processes are kinematically forbidden. As this is not the situation that we assume for the 



present purposes, the corrections that are explained in Ref. [50[ are subdominant in the 
final result for the asymmetry and therefore are neglected here. The results of Ref. 50 
imply however, that a straightforward solution of the Kadanoff-Baym equations should 
indeed regulate the asymmetry when AM 0. Indeed, this is what we explicitly find 
in the present work. 

After these comments, we present the outline of this paper: In Section [2l we use the 
Schwinger-Dyson equations on the CTP, in particular the Kadanoff-Baym equations, 
in order to calculate the asymmetry that emerges from the decays of mixing neutral 
scalar particles. The Kadanoff-Baym equations decompose into constraint and kinetic 
equations. When AM ^ Fd, one may use either of these equations in order to derive 
the asymmetry. If this is not the case, one should use an expansion based on AM <^ 
M. In this limit, the kinetic equations take the form that is familiar from the time- 
evolution of a density matrix in the Hamiltonian approach. In order to gain a qualitative 
understanding, we derive analytical solutions for the mixing propagator in the case of 
two mixing flavours in a constant-temperature background. 



^The Wigner transform is the Fourier transform of a two-point function with respect to the relative 
coordinate, while retaining the average coordinate, see Appendix |^ 
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In Section |3l we apply the same approximation strategies as for scalars to mixing 
Majorana neutrinos. As an additional complication, we need to solve for the various 
spinor components of the singlet neutrino propagator. However, it turns out that the 
flavour dynamics shares the same essential features with the scalar model. Even though 
our main interest is the decay of Majorana fermions, a reader who is less interested 
in these technical details can skip this Section and continue with Section HI There, 
we discuss some features of the analytical solutions for the asymmetry in a model in 
constant temperature background, for both, fermions and scalars. In particular, we see 
that when AM ^ Fd, the effect of oscillations of the singlet neutrinos averages out and 
the standard results for the decay asymmetry from S'-matrix elements apply. However, 
when the condition AM ^ Fd does not apply, the oscillations can have a significant 
effect on the asymmetry. 

In Section |5l we present the modifications of the kinetic equations that arise in an 
expanding Universe. We therefore eventually arrive at an effective theory for Resonant 
Leptogenesis, that is formulated within the CTP framework. In Section El there are a 
short summary of this paper as well as concluding remarks. 



2 Scalar Model 

2.1 CTP Approach to Generating a Charge Asymmetry from 
out-of-Equihbrium Decays 

The perhaps simplest model for generating a charge asymmetry from out-of-equilibrium 



decays is given by |29|, |30|, |46 



The field </:> is a complex scalar that can be associated with a charge, while the Xi real 
singlet scalar fields, that can create a charge asymmetry within when decaying out-of 
equilibrium. We sum over repeated indices z, j. We take here i = 1, 2, for simplicity. The 
generalisation to more than two flavours of x is mostly straightforward, except for the 
analytic formulae involving the diagonalisation of matrices that we present in Section [2^31 
We choose a basis where the mass- matrix M^ is diagonal. The couplings gi are complex 
and of mass dimension one. 

When replacing the gi by Yukawa couplings and the complex scalar by a Higgs and 
a lepton field - which is straightforward to implement, above model describes soft Lep- 



togenesis from sneutrino mixing |64i-|66[| . This may be an interesting phenomenological 
application of the methods developed in the present work. 

In Wigner space, the Kadanoff-Baym equations for the Green functions of the fields 
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Xi and (p are 



A 



e"'^{n^}{A<'>} - e"'^{n<'>}{A^} 



(2) 



= 2e-^^({n>}{A<}-{n<}{A>}) 

and the equations for the retarded and advanced propagators 
1 



A 



R,A 



e-^^{n^ }{ A^'^} ± e-'^im-^}! A«'^} = 1 . (3) 



In these equations, A = A{k,t) and 11 = Il{k,t), and for the particular fields, we 
substitute A — )■ A^^^, U — )■ IIv ,^ a nd M — )■ M^^^. The derivation of these equations is 
discussed in detail in Refs. j48|, |50|-|52|. As a quick reference, we list in Appendix El the 
definitions of the various Green functions and self energies as well as of the o operator, 
that appears in Wigner space. In order to keep the equations compact, we often suppress 
the momentum and time arguments of the Wigner functions in the following 

Next, we keep only the zeroth order contributions in the operator o, which we justify 
in more detail below. By taking the anti-hermitian and hermitian parts, it is useful to 
split the Kadanoff-Baym equations into the constraint. 



4 



A<'> - {M' + n^, A<'>} - {n<'>, A^} = -^[in>, iA<] + ^[in<, iA>] 



and the kinetic equations 

2ik'^dtA<'> - [M\A<'>] 



1 



({in>,iA<}-{in<,iA>}) . 



(4) 



(5) 



Of course, for the single-flavoured field if, the commutator involving vanishes. 
The self-energies for x are given by 

d^k 



and for ip by 



(27r) 



[g,g*iA;\k)iAf{p - k) + g:giiA'^\-k)iA'^\k - p)] (6) 



d'^k 
(2^ 



iAfik)iA^Hp - k) . 



(7) 



We absorb 11^ into a redefinition of M and can therefore effectively drop all terms 
involving it henceforth. 
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Integration of the kinetic equations ([5]) for (^9 yields the evolution of the charge asym- 
metry 



^ip^dtiA<'> = S^ + W^ (8) 



1 f dp° 



2 J 2^ (in>(p)iA<(p)-in<(p)iA>(p)) 

In this form, the equation encompasses both, the washout W^p of the charge of (f as well 
as the source term for the asymmetry S^. When we assume that the charge density of ip 
is small compared to the number density of its quasi-particles, a very useful simplification 
arises (c/. Ref 3^). The washout term W^p can be approximated by including only the 



leading order effects, that do not involve CP-violation. For the CP-violating source S^p, 
we may approximate the field ip as following an equilibrium distribution. Decomposing 

i'^Ax. = iA<f - iA^- (9) 
and using KMS relations fll24p . we then find for the CP- violating source term 

= -9^9- J ^,^S^x, miil) ■ (10) 



where we define 



fl^ = U-^J{g*g,+g,g*). (11) 



The spectral self-energy 11^ is a pivotal ingredient to calculations of the charge asymme- 
try within the CTP approach. In Appendix [HI the concrete expression for its equilibrium 
form is presented. Note that Il-^..{k)/k^ is the decay rate of Xi- 

It remains to determine iSA^.. (q), the deviation of the field x from equilibrium. From 
the constraint equations, we infer that the mass-diagonal components can be written as 

iSA^^Xp) = 2vr5(p^ - Ml) (^(p°)<5/,Jp) + ^(V)5/x.(p)) • (12) 

The functions 6f^ and 6f^ describe the deviation of the quasi-particle distributions of 
the fields x from their equilibrium Bose-Einstein form. Just as the Wigner functions, 
they are time-dependent, but we always suppress the time argument. For the diagonal 
components, the neutrality of the fields Eq- fHOSp . implies that 

SfxM = ^fxuip)- (13) 

Notice that the terms within the constraint equations that involve A^ only contribute 
to the solution for the equilibrium Green function, but not to 16 A^. For a single scalar 



flavour, this matter is explained in detail within Ref. (50 
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It is useful to evaluate the kinetic and the constraint equations in a distributional 
sense, by integrating over some interval of For this purpose, we extract the distribu- 
tion functions from the propagator. The neutrality condition f|T3|) must be generalised 
for the off-diagonal components. In order to keep the notation compact, we therefore 
define 



- { ; (14) 

^fxiP) for P < 

i.e. Sf^{p) depends on p° only through sign(p'^). We emphasise that this is just a con- 
venient definition, because the values of 6f^{p) for far off-shell momenta p are irrelevant 
within our present approximations. We then extract these distribution functions through 

SfW', P) = / ^2/sign(/)i5A^(p) Vp'° I sign(p'°) = ± . (15) 

Here, I± is a small interval around the positive (negative) quasi-particle pole. Note that 
the neutrality condition for scalar fields ( llOSp and hermiticity ( 1106^ imply that 

5/x.(/,P) = ^/x.(-/,p) = 54,(-/,p) • (16) 

Integrating the kinetic equations ([5]) according to Eq. (ITSl) . we obtain 

2ik''dt6f^{k) - [Ml 6f^{k)] = -i{U^{k), 6f^{k)} . (17) 

These equations should be solved by setting k^ to its value at the quasi-particle pole. 
The solution is 6f^(k) for k^ > and 6f^{]<.) for k^ < 0. It can then be completed 
according to Eq. (fT^ . keeping in mind that 5f-^{k) is only physically meaningful for 
momenta close to the quasi-particle poles. 

The off-diagonal components of iA^ can be determined based on two separate ap- 
proximations, depending on the parametric regime. In order to define these regimes, we 
introduce 

rD(p)^max n^^^(p)// (18) 

as a measure of the width of the singlet x o^^? equivalently, its decay rate. Besides, we 
define 

AM=|M^,, -M^,,|, (19a) 

M = ^" ^ — ^ . (19b) 

The two parametric regimes are given by AM ^ Fd and by AM <^ M. When x is 
weakly coupled, as it is a typical requirement in scenarios for baryogenesis from out-of 
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equilibrium decays, we have Fd <^ M, such that there is an overlap of both regimes 
when Tj) <^ AM <^ M. In that situation, both approximations should lead to the same 
predictions for the asymmetry. 

The particular approximation strategies are now as follows: 

• When AM ^ Fd, we can calculate the off-diagonal components of i^A,^ as higher- 
order corrections to the diagonal ones. For simplicity, we assume that only Xi 
deviates from equilibrium. The general case follows from the superposition of the 
procedure that we describe here. In this approximation, the off-diagonal compo- 
nents share the diagonal mass-shells, i.e. the poles of these quasi-particles are 
given by 

i6A^^,ip) = 2n6ip' - MlJ6f^^^ip) , i6A^,M = ^^^ip' - MlJ5f^^M (20) 

for Fd < AM . 

As we show in Section \2.2\ using either the kinetic or the constraint equations, it 
indeed follows that the off-diagonal components of the propagator are suppressed 
by factors of order Fd/AM compared to the diagonal ones. These off-diagonal 
components then enter the source term as given in Eq. (fTOj) . while higher order 
corrections to this are again suppressed by factors ~ Fd/AM. In principle, the 



same approximation strategy is employed in Ref. [30 



When AM <C M, we have no procedure yet in order to determine the precise 
location of the quasi-particle poles, but know that these are located at — M^, up 
to a relative error of order AM/M. We can therefore approximate the off-diagonal 
out-of-equilibrium Wightman function as 

i6A^^^ (p) = 27r5(p2 - M^)6f^^^ (p) for AM < M . (21) 

When using this in order to solve the kinetic equation f[T7|) as described above, we 
expect a relative error of order AM/M. Note that this approximation strategy is 



similar to the one that is employed in Refs. [53|, |54 



In both cases, the result for iA^^^ can then be substituted into the source term f llUp . 

It is now the point to come back to the truncation of higher orders of the o operator, 
a procedure known as the gradient expansion. It is justified, provided we may neglect 
contributions to the finite width of the propagators, which are of relevance when the 



reaction x ^ 2(p is kinematically forbidden (50|. Furthermore, contributions to the 



gradient expansion from fast flavour oscillations, which are systematically treated in 



Refs. [55|-|62[, must not give sizable corrections to the phenomenological results. In the 
resonant regime, such corrections are suppressed by the small oscillation frequency and 
therefore effectively by factors of order AM/M, while in the non-resonant regime, the 
suppression results from the averaging of the oscillatory contribution to the source for 
the asymmetry, as we discuss below. Besides these considerations, the higher order terms 
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in the gradient expansion yield contributions of order H/T ~ T/mp\, where H is the 
Hubble expansion rate, mpi the Planck mass and T the temperature of the Universe, 
which sets the typical energy of a quasi-particle. For temperatures sufficiently below the 
Planck scale, corrections of this type should therefore be negligible. 



2.2 The Regime AM > Td 

In this regime, when using the ansatz fl20|) . we find the following approximate solution 
to the kinetic equations f[T7|) : 

Sf„M=SflSk>-^\ (22a) 

5 fx, (k) = i ^f^^^, ^/x. (k) = -6f^^^ (k) . (22b) 

Xii Xjj 

We assume here that only Xi is excited by deviating from its equilibrium distribution, 
i.e. 5f^ = (^fv L n, (^fv L n = and we set = i^/k^ + . The case with an 

•'Xii •'Xii\t=0^ •'Xjj\t=0 V Xii 

excited Xj can be written down in analogous manner, taking account of the different 
underlying mass shell. As anticipated, the off-diagonal components are suppressed by a 
factor of order Fd/AM. 

When substituting 6f^-. into the Wightman function according to Eq. (!20|) and even- 
tually into the CP-violating source term ( JTOj) for if, we recover the known answer for 
AM ^ Fd and in finite density background: 

S. = -^f 7^2.^(g^ - Aqjig^g* - 9*9,) ^ iJ^imSf^M (23) 

J \ I Xii Xjj 

- Zi[g^ g, g^ ) J (2vr)32^q^ + M'^^^ M^^ - M^^ ^"^^ ' 



cf. Ref. (30 



Taking account of the time derivative f lTTl) and of the anticommutator terms that 
involve the off-diagonal components of we obtain the leading finite-width corrections 
that appear within the denominator: 

ASS AJJ ASS AjJ 

Notice however, that this result still relies on the assumption AM ^ F^. 

The same result can be obtained when using the constraint equations (jlj). For given 
i(5Ajj (and assuming that iSAij = for i ^ j a.t zeroth order in F^/AM), the off-diagonal 
components follow from 

2e - {Ml^ + Ml)] iSA., = m^JSA,^^ + (in^^^ - in^^^) i5A,, . (25) 
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The crucial term on the right hand side arises from the commutator of the spectral 
self-energy and the distribution function and is neglected in earlier literature. Notice 
that the induced off-diagonal correlation is not oscillating. This is in contrast to earlier 



approximations, where the right hand side of the constraint equations is neglected (55 



63|. In conjunction with Eqs. f[T^ and fl2U]) . in particular with the on-shell condition 
/c^ = M^.., we find consistency with the solution ( I22b|) . 

2.3 The Regime AM <C M and Damped 2x2 Flavour Oscilla- 
tions in a Constant Background 

In this regime, the calculation should account for the potential impact of flavour oscilla- 
tions. Scalar and fermionic fields share the same basic aspects of flavour dynamics. For 
the purpose of generalisation, we therefore recast Eq. (|T7|) as 

a,5/ + ^[n,5/] = -i{r,5/}, (26) 

where 

n = Ml/k^ , r = n^/fc° (27) 

and 5f = 6f^. Of course, we can identify this with the evolution equation for a density 



matrix in the Hamiltonian formalism 19, 20, 23-26 



While this is not the case in the Early Universe, it is nonetheless instructive to 
consider the situation when Q and F are time-independent. The solution to equation (l26ll 
is then given by 

Sf{t) = e"^™-5r*5/(t = 0)e5™-|r* . (28) 

Since is diagonal, we have chosen to work in a flavour basis where Q is diagonal 
as well. In this basis, F is non-diagonal in general, which is crucial in order to obtain 
CP-violating effects. Notice that Q and F are both hermitian. Let the matrix 

E = n-iT (29) 

be diagonalised by the transformation 

Sd = . (30) 



We choose the parametrisation (see e.g. Ref. |67|) 




U= , (31) 
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where 



A 

ti 

c 

-Dll 



-22 J 



A+D ' 
1 

Hii + 5, 



D = sign(Re[A])VA^ + Hi2S2i 

^2 - A+D ' 

^ Sl2521 



(32) 



-D22 



A+D ' 
^22 - 5 . 



The matrix 



is diagonahsed by 



= n + ir 



(33) 



(34) 



where V can be constructed in the same way as U. Notice that SJ^ = E^, W = V ^ and 
detU = 1. 

We can thus express 

Sf{t) = U-h-'-°'USf{0)V-'(^-^^'V = U-'e-'-^'U5f{0)U^(^-^^'U^~' . (35) 

Note that Sf{t) is hermitian by construction, provided 5/(0) is. From this formula, we 
can infer the distribution functions of x and their correlations as a function of time in a 
straightforward way. 

For this purpose, appropriate boundary conditions for Sf need to be specified. In 
order to choose these, it is useful to first consider the equilibrium propagators. We 
express these various propagators in terms of their adjugate matrices 



A 



R 



+ in^ v 



adj(A;°^-a;2 + in 



-^adj(/t°^-cc;2-in 



A-^ 

X 



A^ 

X 



1(A^ - A^) 

2^ X ' 



1 



2\V\ 
i 

2\V\ 



adj (A;°' -co^){V + V*) + i adj (U-^) {V - V* 



where 



V = det 



00^ + in 



A 



(36a) 
(36b) 

(36c) 
(36d) 

(37) 



and cu^ = + M^. These solutions can be obtained from the equations for the retarded 



and advanced propagators 



Besides, the equilibrium Wigner functions iA<'^^(/c) 
12 



2A^{k)/{exp{fik°)-l) and iA><='i(A;) = 2A^{k)[l + l/{exi){pk^) -1)] can be consistently 
obtained as solutions to the constraint equations (jl]) as well, cf. Ref. jsoj for a detailed 
discussion. 

As we assume a mass-diagonal basis, where cu^ is diagonal, it becomes clear that 
in the hierarchical limit, this is indeed the appropriate basis for defining the density of 
quasi-particles. In particular, notice that for AM^.^ ^ F, it follows 

iA^.. — 5 ^ . , (38) 

while 

where uf^ = + M^__. Hence, the diagonal contributions are of the Breit-Wigner type, 
as required by the spectral sum rule 



^A:°AJ^ = 5., , (40) 



while the off-diagonal entries are representations of the principal value that cannot be 
associated with particle number. 

Therefore, while it is tempting to choose for Eq. (135|) initial conditions where U6f{0)V~^ 
is a diagonal matrix, because then Sf{t) has no oscillatory contributions, we should in- 
stead choose 5/(0) diagonal in order to describe a state that initially exhibits no mixing 
of the mass eigenstates. For illustrative purposes, we therefore consider a situation where 
Xi is out-of-equilibrium, while X2 is in equilibrium, i.e. 

Eq. (l35l) then yields 

(^(A* + L)*)2ei=mi* - |Fi2pe5-o22*j (^(A + L))^^^^"* - |Fi2pe-^-^22i 



Sfu{t) = 5fo 



5fi2it)=Sfo 



(A + D)^ - IF^-Is 



2 



12 



(42a) 

iFi2(A* + D*) (^e^-Dii* - e^-D22t'^ ({A + Dfe~^2^niit _ |p^^|2g-iHo22t j 



|(A + D)2- |F 



12 



|2| 



(42b) 

5/2i(t) = 5/1*2 (42c) 

|Fl2p| A + Dp (e^^Dll* - e^^D22* I U-^^Dut _ g-fHi522M 

5/22 (t) = 5/0 ■ (42d) 

|(A + D)2- |Fi2|2|' ^ ' 
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In the CP- violating source term for the scalar model f lTOj) . this enters when using Eq. fl2T|l . 

It is interesting to study the behaviour of the solutions fl28|) and fH2|) in the degenerate 
limit, where M^^^ = M^^^^. Because 11^ is real symmetric, is real. We can then see 
from Eq. flTUl) . that the CP- violating source 5';^ is proportional to the imaginary part 
of ^fi2- In Eq- fl^S]) . VL commutes with all terms for AM = 0, since it is proportional 
to the identity matrix. Therefore, no contribution to S^p arises in that situation. From 
Eq. f l42bp . we arrive at the same conclusion, because both A and D are purely imaginary 
for AM = and Re[SDii] = Re[SD22]- Therefore, Sfi2{t) turns out to be real, which 
implies that = 0. 

We should also comment on the case where (A+D)^ — |ri2p = 0. While from Eq. 
it is immediately clear that the result for Sf{t) must be non-singular also in this limit, we 
need to explicitly convince ourselves of the cancellation of singularities in Eqs. fH2|) . From 
Eqs. ([32]), we find that A^ = H12S21 = |ri2p, D = and Hon = Sd22 = Sn - 



in that situation. Substitution into Eqs. (|42|) then explicitly reveals the cancellation. In 
conclusion, we find = for AM = 0, even when the fields xi,2 are distinguishable 
through their interactions gi^2 with the field (f. Of course, one can alternatively arrive at 
the same conclusion through a field redefinition through a unitary transformation of 
that may remove the phases from the couplings g but leaves M^ invariant for AM = 0. 

Finally, for AM ^ Fd, we can cross check the solutions ( 14 2b I) by neglecting the 
oscillatory contributions based on the assumption that these average out in the final 
asymmetry, as we discuss in Section HI In that case, we find 



ind _ _nfi, 



Xii X22 



which is in agreement with the result ( l22l) . where the approximation of non-degeneracy 
has been applied at an earlier stage of the calculation. 



3 Fermions 



3.1 CTP Approach to Leptogenesis 



For our treatment of Leptogenesis, we follow Ref. (33|, where the CTP approach in the 
non-resonant regime is formulated. A lot of that material applies to the resonant regime 
as well, which is the main objective of the present work. Besides, we employ the results 



of Ref. 33 as a benchmark in order to test the results that are derived here. Because 



they are applicable arbitrarily close to the limit of mass-degenerate singlet neutrinos, the 



present results are a generalisation of those presented in Ref. |33[ , and we show that both 
agree sufficiently far away from the mass degeneracy. Notice that we do not reiterate the 
calculation of the asymmetry from the vertex (penguin) correction in the present paper. 

The Lagrangian that gives rise to the standard scenario of Leptogenesis in the un- 
flavoured regime is 



C 
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Here, iVj are the singlet right-handed Majorana neutrinos, that correspond to the fields 
Xi in the scalar model. Their masses are given by the mass matrix M, for which we 
choose a diagonal basis, unless explicitly stated otherwise. Again, we take i = 1,2, for 
simplicity, while a generalisation to more fiavours of N is straightforward. The Standard 
Model Higgs doublet field is 0, and £ is a linear combination of the three left-handed 
lepton doublets of the Standard Model, ipe and ipN are the spinors associated with i 
and A^. Note that the diagonal form of the Yukawa couplings Y, that is assumed in 
the Lagrangian (jH]), can be achieved by a unitary transformation of the left-handed 
lepton fiavours. The SU(2)l indices are contracted in a gauge-invariant way, i.e. (f>ipi = 
—(pA^ABipiB and ip^cl^ = {4'£)a^ab{(P^)b, where e^s is the antisymmetric 2x2 tensor with 
ei2 = 1. The four-component spinor ip^i observes the Majorana constraint (11081) . The 
chiral projectors are Pl,r = (1 =F 75)/2. 

In the Lagrangian of the form (jM]), M may be a general complex symmetric matrix. 
Through field redefinitions of N, it is possible to bring M to a real diagonal form, along 
with an according reparametrisation of 1^. Yet, once a basis is chosen where M is real 
and diagonal in the vacuum, thermal corrections may lead to a real symmetric effective 
mass that adds to M. The off-diagonal components can then be be removed through an 
addtional orthogonal transformation, but in general not in a temperature-independent 
way. Yet, choosing a basis where M is real and symmetric, has the advantage that the 
CP-violating phases are isolated within Fj. Moreover, the separation of the Kadanoff- 
Baym equations into kinetic and constraint equations by taking the hermitian and anti- 
hermitian parts is facilitated when M is purely hermitian. For these reasons, we continue 
the discussion under the assumption that M is real and symmetric. 

The propagators for the left-handed lepton and the Higgs field are listed in Ap- 
pendix O At zeroth order in perturbation theory, the diagonal components of the singlet 
neutrino propagators are parametrised in terms of the distribution functions /tv as 

^S^M = -27r(5(p2 - M^^i^ + Mu) [^(po)/iVn(p) - ^(-Po)(l - fNu{-v))] , (45a) 
'^S^Nu{p) = -27r(5(/ - M^^i^ + M,,) [-^(po)(l - /iVn(p)) + ^(-Po)/iVn(-p)] , (45b) 



(45c) 



iS^M = -^i±MA_ _ 27r5(p2 - Ml){i^ + Mu) [^(po)/7v..(p) + ^(-Po)/7v..(-p)] • 

(45d) 

In general, /at decomposes into equilibrium and out-of-equilibrium contributions as 
In = I'm + ^/a^' where is a Fermi-Dirac distribution. For the equilibrium contribu- 



tions to i^AT, we could then derive finite-width corrections, cf. Eqs. fl36l) and Ref. 50 
for the scalar case. As these corrections are not relevant for the present analysis, we 
do not include them here. Notice that the propagators fj45|) observe the Majorana con- 
straint ffTO . 
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We aim for approximate solutions to the the Kadanoff-Baym equations for the singlet 
neutrinos, which are 



(46) 



, J { ^Nkj } { ^ JV» J { ^Nkj } 



with the self energy 



X 



(27r)4 (27r)4 



(47) 



Through the factor g^^ = 2, we take account of the SU(2)l multiplicity. Since Ref. [331 
omits a discussion of how this self-energy can be obtained from the two-particle- irreducible 
(2PI) effective action approach, that is often used in non-e quil ibrium field theory, we 
briefly present this connection in Appendix [D] (c/. also Ref. [68|). In analogy with the 
discussion of the scalar case, the hermitian part of the self-energy can eventually be 
absorbed in a redefinition of M^, that appears below in the effective equation fl95|l for the 
evolution of the distribution function for the iVj. Technically, this can be performed by 
generalising the methods for Weyl fermions in Ref. 3^ to the case of massive Majorana 
fermions. The contribution that results from tlj^ is real symmetric, but temperature 
dependent. In general, it is therefore not possible to diagonalise the effective mass in 
a temperature-independent way. The temperature-dependent term is given by Eq. (p6l) 
below. 

When replacing the definitions (ITS!) and (fT9|) by 



rD(p) 



^ max 



tr 



(48) 



and 



AM = |Mii - M22I , (49a) 
Ml 1 + M99 

M = " ^ , (49b) 

we can identify the same basic approximation strategies for the calculation of the lepton 
asymmetry from mixing Majorana neutrinos with the calculation of the asymmetry from 
mixing neutral scalars, that is discussed in Section |2l In particular, it is useful to make 
use of the approximation Fd ^ AM, where the crucial off-diagonal components of the 
singlet neutrino Green functions can be obtained as small perturbations to the diagonal 
ones, and of the approximation AM -C M, where the off-diagonal Green functions can be 
obtained from the kinetic equations. Fortunately, since M ^ Fd, due to the smallness of 
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the couplings Y, the parametric regions, where these conditions hold, overlap, such that 
phenomenologically reliable predictions are available throughout the parameter space. 

We should now specify the pivotal quantity pj^. Since for the massless lepton, iS^ 
is proportional to 7-matrices, the same holds true for the self-energy. In particular, we 
may write 

^mjik) = r9^{Y*Y,ti^^P^ + Y,Y;ti^^P^) . (50) 

The expressions for S'^^ r with lepton and Higgs distributions in kinetic equilibrium, i. e. 
for Fermi-Dirac and Bose-Einstein distributions with chemical potentials /i£ and ^(^, are 
given by Eqs. (I114p . When AM ^ is not fulfilled, the off-diagonal correlations of 
the right-handed neutrinos can effectively retain lepton number over an amount of time 
that can be relevant for Leptogenesis. In order to accurately treat the washout of the 
asymmetry, it is therefore important to account for the chemical potentials /i^ and /i^ 
when solving for the Kadanoff-Baym equations for the singlet neutrino propagator 



The extreme case is a setup considered in Ref. [69|, where the right-handed neutrinos 
and their interactions approximately conserve lepton number. On the other hand, when 
AM ^ Fd holds, it is a good approximation to neglect the lepton and Higgs chemical 
potential for the purpose of calculating Pni such that we can set 



All _ y-i^M I s^Afji 

ancH 

Note that |tr[^^'^jj]/A;'' is the decay rate of A^^ 
The one-loop lepton self energy is 



^■^Jk) = ^,ti^g^iY,Y;P^ + F;F,Pl) . (52) 



(2^ (2^^^^)'^'^^ nY:Y,P^iS-J<'>ik')P^iA>'<i-k'') , 

(53) 



which drives the evolution of the lepton asymmetry according to 



S + W = I ^tr \i^>ik)P^iS<{k) - i^fik)P^iS>{k)] , (54) 



d _ „ f d k 

-(n,-n,) 



where ng (ne) is the number density of (anti-) leptons. We decompose this expression into 
the source term S and the washout term W. The source is extracted, when substituting 
equilibrium propagators for iSe and iA^. This approximation is valid, provided the lepton 
charge density is small compared to the number density. In particular, we do not need to 
account for the lepton and Higgs chemical potentials when calculating the source term. 



This definition is related to the similar quantity in Ref. [s^ through = ^^T,'^. 
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The KMS relations then imply, that the source is proportional to the deviation of the 
singlet neutrino from equilibrium, 



S = Y*Y, 



(27r)4 J (27r)4 J {2tt) 



'4 U" 



d'^k 



{2'Kf6\k -k' - k") 



(55) 



X tr [P^i5SM^^{k')P^ (\^<{-k")\S<{k) - \^>{-kyS>{k))] 
d^k' 



(27r) 



rtr 



A 



Pni6SNJk')2Pi^tN{k' 



with equilibrium distributions and vanishing chemical potentials for i and 0, such that 
E'^ is given by Eq. fl52|) . A complete network of equations that determine the evolution 
of the lepton charge in the expanding Universe is presented in Section [51 



3.2 Helicity Block-Diagonal Decomposition of the Singlet Neu- 
trino Propagator 

In order to explicitly demonstrate the consistency of the constraint and the kinetic 
equations when AM ^ Fd and in order to pave the ground for the calculation of the 
asymmetry when AM <^ M, we decompose the fermionic Green functions into vari- 
ous Dirac matrix components. In particular, for a spatially homogeneous and isotropic 
problem, it is useful to notice that the helicity operator h = A;*7°7*7^ commutes with 
the various Dirac matrices that appear in the Kadanoff-Baym equations (HSI), such that 



helicity is a good quantum number. This suggests the decomposition [51 

1 



71 



i6S 



N 



h=± 



Nh 



Nh 



-(l + /iA;V)®p%;, 



(56) 



where h = ± denotes helicity and a*, p* are Pauli matrices. Note that the functions Qah 
are understood to be hermitian matrices in flavour space. It may be useful to compare 
with Eq. f ll29p . in order to acquire a quick understanding of the functions Qah- 

Multiphcation of the Kadanoff-Baym equations fH6|) by {1, —^17*7^, —ihY, —7^}, tak- 
ing the Dirac trace and truncating the diamond operators yields the equations 

/ 9oh\ 



[j^h-Y*Y'ChL-YY^ChR\ 



9ih 

92h 

\93hJ 







(57) 



where Y* = (^1,^2), 



( 



ik" + \dt 


\M 





i/i|k| 


\M 


-ik° + \dt 


/i|k| 








-^|k| 


-iko + \dt 


M 


i/i|k| 





-M 


-ikO + 



(58) 
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and 

C,L,R=yX (59) 
~^JVL,R -F ^^*^'iVL,R -FiSjvL,R ~ i^^*^'5vL,R 

±iE;;^L,R + i^^'^L,R ~^jvl,r T ^^*^jvl,r 

\ ='=^iVL,R + ^^'^iVL,R ^^iVL.R T ^^*^JVL,R / 

The kinetic equations correspond to the hermitian part of Eqs. ([57]), which is 

+ i[M, ^71,] = 5^ -y (S^°L,R ± /^^^^^l^) {Tl,r, g^h T Qzh} , (60a) 

g^n + 2/i|k|(72/. + i[M, g,^] = y (e^°l,r ± /^^'^^l,r) (-{Tl,r, Qih} T i[TL,R, ^72^J) , 

(60b) 

g^u - 2/i|k|(7i, + {M,g^u} = y (s^°l,r ± hUt^^^^ (±i[TL,R, - {T^,K.92h}) , 

(60c) 

9ZK - {M, = E - Y (^^L,R ± ^^'^^Ir) {T^l.r, 9zh T 9oh} , (60d) 



± 

where we define the shorthand notations 



Tl = Y*Y* , Tr = YY^ . (61) 
The constraint equations are the anti-hermitian part of Eqs. flFTI) and read 

-2iA;°^7oh, + 2ih\k\g;h + i{M, g^h] = -y {^il,K ± /^^^^^I^^r) [Tl,r, (7oh T 9^h] , 

(62a) 

-2iA;°<7i, + i{M, = X] y (s^°l,r ± hUt^^^^^ (62b) 

X (-[Tl,r,5'i/.] =F i{TL,R,fi'2/.}) , 
-2\k'g,H + [M, <73.] = 5^ y (s^^r ± /^^'S^I,r) (62c) 

X (±i{TL,R,fi'lh} - [TL,R,5'2h]) , 

-2ifc°(73h + 2i/i|k|(7o/. - [M, g2H] = -y (s^°l,r ± /^^^S^I^^r) [Tl,r, (^sh T ^?oh] • 

(62d) 
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3.3 Derivation of the Asymmetry for AM ^ Fd 

In the regime AM ^ Fd, we calculate the off-diagonal components of iSS^ij as small per- 
turbations for a given diagonal out-of-equilibrium distribution. We verify that constraint 
and kinetic equations give consistent answers. This is in analogy with the procedure that 
is presented in Section [2l2] for the scalar case, and the error is controlled by the expansion 
parameter Fd/AM. 

We first derive the source of the lepton asymmetry from the constraint equations ( 162|) . 
Assuming a mass-diagonal basis and neglecting the collision term, we obtain from the 
diagonal components of Eqs. (16 2 p the relations 

9ihii{k) = ^gohii{k) , g2hii{k) = , g3hii{k) = ^^gohuik) . (63) 

For these diagonal components, we assume helicity-independence, i.e. go+u = go-u- 
This follows when assuming that there is no initial helicity asymmetry within the singlet 
neutrinos and no initial charge asymmetry within the lepton and Higgs sector (and hence, 
no CP-asymmetry). Because the tree- level decays and inverse decays then produce 
no helicity asymmetry, also the non-equilibrium distributions for the Ni are helicity 
symmetric at leading order. A helicity asymmetry in iVj is only incurred through CP- 
violating loop effects and the back-reaction of the produced lepton asymmetry, which 
results in sub- leading contributions to the non-equilibrium distributions for the N^. Note 
that symmetric initial conditons are naturally established within the parametric range 
of strong washout. Through Eqs. we may relate gohu to the distribution functions 
as 

sign(A;°)2y/k2 + Mf^6f^M27r6{k' - M^) = g.^^k) . (64) 

From the expression for the source term ( 155|) . we see that we need to determine the off- 
diagonal components (that is i 7^ j) of P^i6SNijPi^. Therefore, we employ the diagonal 
solutions to derive the off-diagonal contributions as perturbations of order F^/AM. 
Because AM ^ F^, lepton number is rapidly violated once the inverse decay into a right- 
handed neutrino has occurred. We can therefore neglect the lepton and Higgs chemical 
potentials when calculating the singlet propagator and use the approximation (ISTl) . The 
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constraint equations for the off-diagonal components {i ^ j) are then given by 
-2ifc°^o^5 + 2ih\k\gl,l + i{Mu + M,,)g\S, = \ {^^H ± hUT.i\^i {g,^u T 9,^^ , (65a) 
-2ik'g\i]^^ + i(M,, + M,,)glj]^^ = 1 {T.i% ± hUT.i\^ {g^^^ t i92hu) , 

(65b) 
(65c) 

-2ifc°53^5 + - - M,,)gi,l = tI {^il ± hk^Ei]^) {go,u T 9,hu) , 

(65d) 

where g^^^^^i^ = Y*Yjg^,^^ij, gf^^^^^ = Y*Yigf;^^ij and go^stij = 9o,3hij + 9f,3hij- Because of 
hermiticity, fifo shji ~ (fi'o'shij)*- assumed here, that only ghu ^ 0, while ghjj = 0. 

The general solution can be constructed by superposition. 

Given the diagonal solutions fl63|l . we can now straightforwardly determine the gaUj 
for z 7^ j as solutions to the system of linear equations fl65|) . In particular, the flavour 
off-diagonal vector- and pseudovector-densities are 

9o,u = fcO(M,^, - Ml,) [^^ ^^^^^) ' ^^^"^ 

We have now assumed that only A'^i deviates from equilibrium. As stated above, the 
general case, where A^", deviates as well, can simply be obtained by superposition. 

Note that these sub-leading off-diagonal solutions do not obey the leading-order sym- 
metry property, according to which g^h is helicity-even whereas g^h is helicity-odd. For 
the CP-violating source, above results combine to a compact form that observes the 
symmetries of the problem: 

y:y,p^i5Sn,, {k)PL = --y:y,\ [p^^'p^ - ^^r^^t^^l] {9^^, + 93M,) (e?) 

h=± 

--[Y, Y, -Y,Y, )p^2^_^2/R7 Pl^^n,9ohii 

= -{YfYi - Y^Yf)i^^^^pXP^'^5fm.W-5{k- - Mf) 

Since 7V2 is assumed to be in equilibrium, 5fN22 = 0. Notice that only contributions from 
^o/iii enter here, while those from cancel since they are multiplied by the conjugate 
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Yukawa couplings. Substituting Eq. f l67|l into Eg. fl55l) . we obtain the known result for 



the CP- violating source of the lepton asjTumetry [33 



which includes the finite-density corrections. We refer to this expression as the perturba- 
tive result for the source term, as it can be calculated as a perturbation to the diagonal 
singlet neutrino propagator or, alternatively, from a loop expansion of Feynman diagrams 



as in Ref. 33 



While above result is obtained from the constraint equations fl62|) . we now check the 
consistency with the kinetic equations (1601) when using the same approximations: When 
AM ^ Fd, the diagonal (ll)-densities track the deviation of A^i from equilibrium, which 
provides the thermodynamical breakdown of time-reversal invariance, which is necessary 
for baryogenesis. Compared to these diagonal components, the off-diagonal ones are 
suppressed by a factor of order Fd / AM. In the absence of fiavour oscillations, the time- 
derivatives of the off-diagonals are additionally suppressed by the same parameter, such 
that we may neglect them. With these approximations, the resulting set of equations 
obtained from the (12)-components of Eqs. (HI]) is 



i(Mn - M22)^l;^2 = ± hUT.i\^){gohii T Qshii) , (69a) 

2h\k\glj^'^, + i(Mn - M,,)g^/^^, = -^(S^la ± hk^J:i\,)ig,j,u T mhu) , (69b) 

-2h\k\g\S, + (Mn + M22)glS2 = T^K^i^ ± hk^J:-^,,){g„,n T mhii) , (69c) 

-(Mn + M22)^72';?2 = ±^(S^i2 ± hyj:i\,){gohu T 93hii) • (69d) 
Indeed, from these equations, we can straightforwardly reproduce the result fl66|) . 



3.4 Asymmetry for AM <C M 

In order to obtain a solution for the CP-violating source also in the regime where AM 
is of the same order as Fd or smaller, we could aim for a full solution of the kinetic equa- 
tions f l60|) . Within an intuitive picture of flavour oscillations, the meaning of Eqs. fl6Ql) 
is yet somewhat obscure because of the mixing of the gah, which is a result of the spinor 
structure. A diagonalisation of this spinor structure is not possible in general, because 
the dispersive contributions from M and the absorptive ones from cannot be diago- 
nalised simultaneously. 

In the close-to-degenerate regime, where flavour oscillations are important, crucial 
simpliflcations can however be achieved. We do not need to specify the locations of 
the quasi-particle poles of the diagonal and off-diagonal components of the propagators 
exactly, but only make use of the information that these are located at = M^, up to 
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a relative error of order I\MjM. We can then integrate the constraint and the kinetic 
equations fl60|) over an interval X of length of order AM over the region of the quasi- 
particle poles. As a result, we obtain again the kinetic equations f l60|) and the constraint 
equations (1^ with the replacement gahij — ^ ^fahij, where 

• . On f ^ i Sfah^j{p) for + 
— sign(p )gahij[p) = \ f r ' (70) 

^TT 5/afeij(p) for - 

and Xj- is again an interval around the positive or negative energy quasi-particle pole, 
respectively. The error results from the averaging over values of in the constraint equa- 
tions and can be estimated to be of relative order AM/ M. In order to keep the notation 
compact when we have to distinguish the cases of positive and negative frequencies, we 
adapt the pragmatic definition (fT4ll to the fermionic case, 

rr ^ . j Sfah^j{p) for / > 

[ ofahijip) for p'' <0 
When Td < M and AM < M, it follows from Eq. fITO or from Eqs. §^ that 

SflH^,ik) = 5^,,,ik)^^^g^ (72) 

and from Eq. (162bp that 



ShH^,{k) = 6foM,{k)^^^^±^. (73) 
Using these approximate constraints, Eq. (l6Ual) becomes 

SfoH + ^i[M^ Sfon] = - E Y {^nIk ± /^^^S^Ir) (l T ^) {Tl,r, Sfon} (74) 



- ^7,.<| Re[Y*Y*]^ - iMm[Y*Y*]^, 6fo, 

where k = (|k|, k^k) and in the last row we have neglected the lepton and Higgs chemical 
potentials according to Eq. f lSTj) . As discussed above, the latter approximation is not in 
general justified when AM < Fd. The off-diagonal components of the out-of-equilibrium 
Wightman function for the singlet neutrinos can then be approximated as 

g^Mjip) = 27r5(p2 - M^)2p''6faMAp) > (75) 

which is accurate up to a relative error of order AM/M. In order to calculate the 
CP-asymmetric source, this should then be substituted in Eq. ( !56l) and eventually in 
Eq. (EH]). 
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We can easily perform the consistency check that Eq. fl7^ reproduces above results 
when Td < M and AM < M but AM > Td- Neglecting the lepton and Higgs 
chemical potentials according to Eq. f l5T]) . the off-diagonal terms that arise from the 
out-of-equilibrium singlet neutrino A^^i are 



0hl2 



2i 



3/il2 



2i 



Ml -Ml 2 '^'^ 

k^ + h\k\ i^^Q 
Ml - Ml 2 

k^-/i|k|A:° g^ 
Ml) 2 
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(76a) 



5/. 



O/ill 



A;0(M2 
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hk'T. 



(76b) 



k2 + /i|k|A;0 



11 



MI2) 2 



5/( 



O/ill 



where we have neglected the time- derivatives on the left-hand side of Eq. (|74l) . Therefore, 
these equations cannot take into account possible flavour oscillations for the iVj, which 
should however be irrelevant for the lepton asymmetry when AM ^ Fd, as we explain 
in Section |H Note that these expressions can as well be obtained from Eqs. fl66|) when 
taking Mn — )■ M22 and using the definition flTOj) . As a further consistency check, we also 
note that in the regime where AM ^ Fd but AM <^ M, we can exchange Mn -H- M22 
within Eqs. (ES]) and Mn + M22 O 2Mii,22 in Eqs. ([60]). Within the result ([66]) for 
(70/112 and g3hi2, this indeed only incurs a relative error of order AM/M, which is a small 
correction. 

Now, when we relax the requirement that AM ^ Fd, Eqs. (!74l) can be straightfor- 
wardly solved numerically. Besides, it is also instructive to obtain analytical solutions in 
a situation without expansion of the Universe. These can be obtained when substituting 



n = — 



M^ 



and 



yAO _|_ uv.ty. 



'^Ai 



It 



h\k\ 



(77) 



h ■ y-^ 



h ■ y^ 

ihlm[Y*Y']^-=^ 

k^ 



in Eq. f[2H]) . In the last row of Eq. f[77)) we have again neglected the lepton and 
Higgs chemical potentials according to Eq. ( 151]) . Before we discuss some features of 
these solutions in Section HJ we note an additional consistency check. When neglect- 
ing oscillatory contributions, as it is appropriate for Fd <^ AM, we can immediately 
recover the perturbative solution fl76ap from Eq. (142bl) . Eq. (176bl) is then obtained 
through g^h = (^|k|/^°)fi'o/n which also holds for the flavour off-diagonals in the close-to- 
degenerate regime, as we have demonstrated above in Eqs. (I72|73l) . When AM ^ Fd, 
we therefore consistently reproduce the result for the lepton number violating contri- 
bution to the off-diagonal singlet-neutrino Wightman function f[67]) and eventually the 
well-known result for the CP- violating source 
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We have therefore estabhshed that for AM ^ Fd, the kinetic equations f l26|) accu- 
rately reproduce the classic result for the resonantly enhanced lepton asymmetry [s], [g]- 
10|, that arises as the zero-temperature limit of Eq. fl68|) . In order to achieve this 
agreement, the helicity structure exhibited in Eqs. f in|) and f l77|) is crucial. For the 
lepton asymmetry from decaying Majorana neutrinos, to our knowledge, this is the first 
time an explicit connection between results from the S'-matrix approach and from kinetic 
evolution equations is established. 

4 Illustrative Features of Analytical Solutions when 
Ignoring Washout and the Expanding Background 

We specialise here to the regime where the approximation AM ^ M is valid and first 
consider the scalar model. Therefore, we use the form (1211) for the off-diagonal Wightman 
functions of x- The source term for the asymmetry flTUl) can then be expressed as 

When adding the positive and negative energy contributions, it is useful to notice the 
transformation properties of the quantities in Eqs. 



A — > -A*, (79a) 
D — > -D\ (79b) 
Soii ^ — > . (79c) 

As a consequence, 5f^^^{k^,'k) = 6f^.-{—k^,'k), a relation that can also be derived from 
the neutrality condition (llOSp imposed on iA^. Use of the result f H2|) then yields 



^ = -i{gfg;' - gl'g'^) / (5/o(q)- ^ ^ (80) 

' (27r)32v/^^TM^ '°'^'|(A + D)2-|n^^^/a;(q)|2|^ 



X 



(A* + D*) (^e5=Dn* _ e5=D22*j (^(A + Dfe~-^-^^'' - \I[-^Ju{ci)\^e--^^^^^'^ + c.c. 



where we have used that n^(A;°, k) = — n^(— k), w(q) = \/ + fvP and we evaluate 
= n"^(ci;(q), q). The parameters A, D and H are given by Eqs. (I27|32p . 
In order to assess the role of the oscillatory contributions within this source term, we 
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integrate 



oo 



X 



\{A + Dy-\U^Ju;{ciW 



(81) 



|A + Dp(A + L)) 



roil |(^Dii — ^022) — i(rDii — r 



- (A* + D* 



X12 



2V-iJii - D22, 
1 



H^Dll — ^022) — 2(-'^Dll — rD22) rD22 



C.C. 



where we define Qdu = RefHon] and Tdh = — Im[HDij]. This result corresponds to the 
asymmetry from inverse decays while Xi approaches equilibrium, when washout of the 
charge in ip is ignored. When AM ^ Fd, this becomes 



-2itef 



*2 2\ 

9i 92) 



(27r)32Vq2 + M2 



5/0 (q) 



a;(q) 



M2 - M2 n-^ 

Xll X22 -^-^Xll 



(82) 



The same answer can be obtained from integrating the source (!23|) for ^ AM, when 
substituting 5/^^. = (5/oexp(— (n-^j^(ci;(q), q)/a;(q))t). We can therefore conclude that 
indeed, oscillatory contributions to the CP-violating source approximately cancel in this 
regime. 

Next, we compute the source fl55l) for the lepton asymmetry in a static background. 
For this purpose, we need to calculate the relevant components of the singlet neutrino 
propagator, P-Bi5SNij- These can be expressed through the distribution functions ( l42|) . 
where the parameters A, Z), F follow from Eqs. fl52]) and (1771) . As we ignore the effect 
of washout, we also neglect the lepton and Higgs chemical potentials when computing 
the singlet propagators and make use of the approximation ( 15T|) . Notice that there are 
now two separate solutions for the helicity states /i = ±, and we can therefore identify 
the solutions ( l42l) with 5/o/i, ^-e. 5fohij{k) = 6fij{k). The solutions thus following 
have the property Sfohij{k^,\i) = 5/o/jjj(— k) = 6fohji{~k^,^), as a consequence of 
the Majorana condition (11091) and the hermiticity (II 101) of the distribution functions. 



Making use of this relation, we can conveniently express the sum over positive and 
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negative frequencies that enters into the source for the asymmetry fl55l) as 



5(k) 



Y*Y 



2tx 



-tr 



,A 



(83) 



tr 



fcO=±w(k) 

h=± 



sign(fc°)5/, 



Ohij 



h\k\ 



h=± 



^|k| 

h\k\ 
~k^ 



{5 f Ohij + Sfohij) 



Qhij 



Ohij ) 



fcO=a;(k) 



/i=± 



A; ■ Ysjq 



{5 f Ohij - 5 f Ohij) 



k ■ Sjy 



k^ 



{Sfohij + Sfohij) 



fc"=i^(k) 



where we have used that t^%k^,k) = t^%-k^,k) and t^'{k^,k) = -tf{-k'^,k). 
Above result for S(k) should be useful when 6f is calculated numerically, i.e. for phe- 
nomenological studies of Leptogenesis in the Early Universe, as we outline in Section |5l 
When we neglect the expansion, we can again use the analytic solutions for damped 
flavour oscillations in Eqs. f l42p . Using Eqs. f l77p and f l32p and defining 

(84) 



Sfl2 = 5/i2/(iri2) , 
we obtain the relations 

5/+oi2(c^(k),k) = -E^2ih = +, A;" = a;(k))5/;2(a;(k), k) 

5/„oi2(u;(k),k) = El,{h = +,k' = a;(k))5A2(a;(k),k) , 

5/+oi2(-u;(k),k) = -El,ih = +, = uik))5f*,iuik),k) 

(5/_oi2(-u;(k),k) = E,2{h = +,k' = u{k))6f*,iu{k),k) . 



(85a) 
(85b) 
(85c) 
(85d) 



Note that in the mass-diagonal basis, Eij = —iTij, for i ^ j. Above factorisation can 
easily be understood when noting that |ri2p, as given by Eq. (1771) . does not depend 
on helicity. Moreover, from Eqs. ( l32l) and (1771) . it follows that Ed does not depend on 
the helicity either. Therefore, the complete helicity dependence is isolated within the 
explicit front factor of ri2 in Eq. fl42bp . Substituting Eqs. fl85|) into Eq. fl83l) . we obtain 
the simple expressioiJl 



S{k) = i{YrY^ - n'rr)2^7.^^^^E-^E7(5A2 + 5/1*2) • (86) 

by the sum convention is 



■^Note that the summation over i,j = 1,2 that is understood in Eq. 
exphcitly performed in the following equation. 
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Figure 1: The source iS(k) for the lepton asymmetry over time t for Mi = 1, M2 = 1.1, 
T = 1, Yi = 0.1, I2 = 0.2 + O.li, i.e. AM 3> Fd, and initially a vanishing distribution 



for A'^i and an equilibrium distribution for N2. The solid line is the result (!83|) . the 
dashed line the standard perturbative limit ( 187|) . In the left panel, we take |k| = 0.5, 
and consider the asymmetry resulting from a non-relativistic singlet neutrino. In the 
right panel, the singlet neutrino is relativistic, |k| =5. 
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Figure 2: The source S(k.) over time t for the lepton asymmetry for Mi = 1, M2 = 1.01, 
T = 1, Yi = 0.1, Y2 = 0.2 + O.li, i.e. AM ~ Fd, and initially a vanishing distribution 
for A'^i and an equilibrium distribution for The solid line is the result (!83|) . the 
dashed line the standard perturbative limit ( 18 7p . In the left panel, we take |k| = 0.5, 
and consider the asymmetry resulting from a non-relativistic singlet neutrino. In the 
right panel, the singlet neutrino is relativistic, |k| =5. 

This is the source for the lepton asymmetry that results from inverse decays of A'^i while 
it approaches the equilibrium distribution in a background of constant temperature. As 
a consistency check, the known result for AM <^ M but AM ^ F^ is recovered for 
^/i2 = ^/o[i^(k)/(M^^ — MI2)] exp(— Flit), which results from neglecting the oscillatory 
contributions in Eq. (142bp . 

We now compare the result from Eqs. fl83|) or fl86l) to the standard perturbative 
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approximation 



2V^TW Mf, - Ml/-^^'^^^ ' ^^^^ 

For this purpose, we use the solutions (l42l) with (5/o(k) = — /'^^(k), which corresponds to 
a vanishing number of A'^i at t = in the particular mode k, while N2 is in equilibrium. 
The other parameters are given in the Figure captions. Their choice does not correspond 
to a systematic study of parameter space, but the examples should be illustrative and 
representative. In Figured! we illustrate the situation when AM ^ Fd. The oscillation 
frequency is much larger than the damping rate, such that we may expect that these 
oscillations average out in the final result for the asymmetry, see the discussion above. 
We choose examples of non-relativistic and relativistic singlet neutrinos, because in the 
relativistic regime, the result turns out to be the remainder of an incomplete cancellation 
of larger terms involving E;^° and E'^', which is therefore explicitly verified. The results 
in Figure [2] correspond to a point in parameter space where AM ~ Fd- Since the 
damping rate is now of the same order as the frequency of flavour oscillations between 
A'^i and these should not be neglected when calculating the asymmetry. 

This brings us to comment on earlier expressions on the behaviour of the asymmetry 
when approaching the resonant regime. This matter is discussed in Refs. 0, [sl, [q], |47| . 
where expressions for the asymmetries are given. Oscillations are neglected there, an 
approximation that is not always suitable, cf. Figure |2l It is however emphasised 
in Refs. 0, S, H Q, that in the regimes AM ~ Fd and AM < Fd, the standard 



perturbative result breaks down and a suitable resummation technique for the singlet 
neutrino propagator should be employed, which is essentially what is performed within 
the present work. It has already been mentioned in Section 12.31 that when choosing 
U6f{0)W as a diagonal matrix, there are no flavour oscillations. It should however be 
emphasised that this corresponds to a peculiar initial state that already bears a CP- 
asymmetry. From Eqs. f l32p and (15^ it then follows that 

5(k) = i{YfYi - Y^Yf)2g^,^^^ti^ti^ + 6fo,u • (88) 



Notice that A + D 7^ 0, as long as Yi 7^ 0. This result is in agreement with Refs. [45|, |46 



where the time evolution of mixing scalars in the vacuum is considered. It generalises the 



findings of Refs. [45|, |46| to mixing Majorana fermions, that decay in a finite-temperature 



background. In the limit AM S> Fd, the term ~ S12S21 in the discriminant of D can 



be neglected, and the results of Refs. [4, |47[ for Majorana neutrino decay in a zero 



temperature background is recovered. We emphasise however, that when the condition 
AM ^ Fd does not apply, the full time-evolution of the off-diagonal components of the 
singlet neutrino propagator has to be calculated, because due to oscillation effects, these 
are not simply proportional to the diagonal components in general. 

Finally, let us briefly discuss the limit IF12I -C iFii — F22I and IF12I < AM in order to 



facilitate the comparison with the results in Ref. (72|. In this limit the term ~ S12S21 in 
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the discriminant of D can again be neglected, and Eq. fl86l) for the CP- violating source 
reduces to 



5(k) 



MiiM22{M^i - Mr 



22) 



X 



-Flit 



cos(Aci;t) 



Tii-r 



22 



2 Aw 



sin(Aa;t))e-(^"+^-)*/2 



5fo 



(89) 



where Aw(k) = \Mf^ - M|2|/(2w(kyh This result is found to be in agreement with the 
analytical results presented in Ref . j72| Furthermore, the structure of the denominator 
{M^i — M22Y + i^^(rii — T22Y agrees with the results in ref. 47|]. However, we want to 
emphasize that the parametric regime with |ri2| ^ iFn — r22| is by no means generic in 
resonant leptogenesis. 



5 Effective Theory in the Expanding Universe 

In this Section, we present a network of equations that describes the simplest version of 
resonant Leptogenesis in the expanding Universe. Equations that are valid sufficiently 
far away from the resonant limit, i.e. for AM ^ Fo, are presented in Ref. 33|. In the 
present context, we need to generalise Eq. ( FMl) in such a way, that it encompasses the 
effect of the expanding background, in particular how this induces a deviation of the 
distribution functions of the singlet neutrinos from equilibrium. 

First, we transform the time coordinate t to the conformal time 77 through the re- 
lation dt = a{rj)dri, where 0(77) denotes the scale factor. In the radiation dominated 
Universe, 0(77) = aurj, where cr, is a constant. The scale factor is related to the physical 
temperature by 

1/4 



OR'Tipi / 45 



where g^, denotes the number of relativistic degrees of freedom. Note that the comoving 
temperature Tcom = o,{v)T is constant. There remains a freedom of parametrisation for 
a-R. For kinetic equations in the early Universe, z = Mrj is a convenient choice, where 

M , , 

^ = ^, (91) 

and therefore 



2 V 9* 



"'Note that in Ref. [7^ the initial state corresponds to vanishing number densities for both neutrino 
flavours A^i and 7V2. Furthermore, the analytical results in (tH assume the non- relativistic limit Mi ^ T, 
such that Sj^^(fc) « fc''/(327r). 
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The choice of the conformal time t] instead of the comoving proper time t has the 
advantage that the Kadanoff Baym equations are of a form that is very similar to the 
one they take in Minkowski background. The only modification is that the mass terms 
need to be rescaled as M — )■ a{rj)M and the momenta are to be understood as comoving 
momenta. The derivative with respect to 77 [or, equivalently, z/M^ when using Eq. fl^?2]) ] 
is denoted by a prime. 

Using these parametrisations, one can obtain a term that represents the deviation of 
the singlet neutrino from equilibrium, that is induced by the expansion of the Universe. 
When AM <^ M, the equilibrium distribution for the singlet neutrinos is 

and f^ijiz) = for i ^ j. When we are not working in the mass eigenbasis, this 
relation is only approximately correct and there will be off-diagonal elements that are 
suppressed by a factor of AM/M when compared to the diagonal ones. While this is a 
tolerable error, it could easily be corrected for by diagonalising M through an orthogonal 
transformation. The expanding background induces a decrease of the temperature, such 
that consequently. 



and fN-j{z) = for i ^ j. 

The generalisation of Eq. ( |74l) . valid in Minkowski background, to an expanding, 
radiation dominated Universe is now given by 

(95) 



where = ± ^J\^^ + a{7]fM^ . 



Note that this equation makes no assumption about the form of M, except that it is a real 
symmetric matrix, which allows to account for temperature-dependent effective mass- 
corrections that may be incurred through Eq. (l95l) is therefore flavour- covariant 
with respect to orthogonal transformations of Ni. Allowing for a complex symmetric M 
would lead to a more complicated form of Eq. fl95|) . but it would not capture addtional 
physical parameter space, as the complex phases in M can always be absorbed in Y , cf. 
the discussion in Section 13.11 Therefore, Eq. fl95|l can be used also in situations where 
the mass eigenbasis is time-dependent due to finite-temperature effects. This should be 



31 



of importance close to the resonant limit, since the thermal corrections to M that arise 
from are given by (for momenta |k| ^ jM^^'^l, which covers the relevant contributions 
of phase-space) 

Mj^' = ] {y;Y, + Y^;) , (96) 



which is of the same order as k^Tj). 

Assuming that the lepton charge density is small compared to the number density 
(or, equivalently, the lepton chemical potential fie is small compared to the temperature 
T), we extract from Eq. fl5^ the washout term as 



d k f d k f d k \4 r-4 / 



X 



tr [(PRi5>,,.(A;')PLiA<(-A;") - PRi5<,,.(fc')PLiA>(-A;")) i55,(A;)] . 



The deviation of the lepton propagator ibSnilz) from the form with vanishing chemical 
potential can be expressed through 

Sel'^l/'^ el'^l/^ 
5/,(k) = -5/,(k) = (n, - ne) X ^3 ^ = , (98) 

where 

\8Se,(k) = -PL|^27r(5(fc2)(5/f (k)sign(fc°) . (99) 

In contrast to the washout term quoted in Ref. jssj, we include here the off-diagonal 
elements of the singlet neutrino propagator. When AM ^ Fd, these are suppressed 
compared to the out-of-equilibrium components by a factor of Fd / AM, i. e. we can safely 
neglect these in the non-resonant regime. When AM is of order Fd or smaller, neglecting 
the off-diagonal correlations is not necessarily a good approximation. Especially in the 
weak washout regime 6foh may be of similar size as f^^, and these corrections may be of 
importance. In the strong washout regime however, the deviation of the singlet neutrino 
density from equilibrium is small, such that we can substitute the equilibrium propagator 
for iSn in Eq. (197|) . which is diagonal in the mass eigenbasis, when neglecting finite- width 
contributions. 

The collision term on the right-hand-side of Eq. (p5|) is valid when the chemical 
potentials of the Higgs and lepton fields can be neglected. Once this is not the case, it 
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should be replaced with 



S>° + hk^T>:^) [l - ^) { Y*Y*, /o4 (100) 
2 ^^>A - hk^^>l^) (l + ffl) {YYt, 
+ A (S<° + hUi:<l) [l - ^) {Y*Y^ 1 + /o4 

+ \ (S^°R - /^^^S<g + ^) {YYt, 1 + . 

Here, /oh = /g^j' + 5/o/i is the full distribution matrix of N , that also includes the equilib- 
rium distribution /^"^(k) = diag(l/(exp(v/k2Til^/T) + l),l/(exp(v/k^TM^/T) + l)). 
The S^'^ can be obtained from the expressions for Sj^ in Appendix [B] when using the 
Kubo-Martin-Schwinger (KMS) relations 

SkR(^) = -e('="^'^'^^'^^)/^S>L,JA:) . (101) 

It should be useful to expand this term in the small ratios fie^^/T, in order to separate it 
into a contribution of the form as in Eq. (l95l) and an additional one, that only depends 
on the chemical potentials. A detailed study of this matter is subject of ongoing work. 

Of course, Eq. ( 15^ must be transformed from comoving to conformal time as well, 
with the result 

-^ini-ni) = W + S. (102) 
ari 

For AM ^ Fd, the procedure of solving for the lepton asymmetry is described in detail 
in Ref. 33|. We first determine Sfohij{k) from Eq. fl95l) and the approximation fl5T|) of 



vanishing chemical potentials for Higgs and leptons. The result is then to be substituted 
into the source term fl55l) and the washout term (197|) . When AM < F^, we should not 
expect the approximation flM]) substituted into the kinetic equations for the oscillating 
singlet neutrinos ([25]) to be sufficient, as it may fail to accurately describe the washout 
of the lepton asymmetry. In that situation, the evolution of the lepton asymmetry (11021) 
and of the correlation of the singlet neutrinos (!95|) should be solved as a coupled set 
of equations. In principle, the lepton and the Higgs chemical potentials are further 
processed by spectator effects. A good approximation for first phenomenological studies 
may be to consider an unfiavoured regime, where the charge in £ is conserved, up to the 
interactions with iVj, that are accounted for already by our equations. The hj^ercharge 
asymmetry within (p should get distributed over a large number of remaining Standard 
Model degrees of freedom, such that we may set fi^ ^ within a first approximation. 

Note that for the source term, the form given by Eq. fl83l) may be particularly useful, 
because we only have to calculate the positive frequency solutions for Sfohijik), while 
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the negative frequency solution follows as a consequence of the Majorana condition. 
Eventually, Eqs. f l95|) and (11021) can be integrated in order to obtain the final lepton 
asymmetry. A numerical study of the evolution of the lepton asymmetry, using the 
methods that are outlined in this Section, should be performed in future work. In 
particular, when AM ^ Fd, it appears plausible that for an equilibrium deviation 
induced by the expanding background, the oscillating contributions visible in Figure [1] 
vanish altogether, because the deviation of the singlet neutrinos from equilibrium is 
generated continuously throughout Leptogenesis. Deviations that appear at different 
times oscillate and have different phases, such that the oscillations may effectively vanish 
due to decoherence. However, when AM ^ F^ or AM <^ F^, the methods outlined in 
this Section may be suitable to make predictions for resonant Leptogenesis that are not 
accurately covered by the perturbative result f l68|) for the source term. 



6 Summary and Conclusions 

In this paper, we have developed the theory of Leptogenesis from mixing Majorana 
fermions and from mixing scalars in the CTP formalism. Our calculations rely on the 
expansion parameters Fd <C AM or AM <C M. Due to the smallness of the couplings, 
Fd ^ M, such that the parametric regimes where the two approximations apply overlap. 
When Fd <^ AM, the usual result for the CP-asymmetry from the wave-function cor- 
rection applies js], IgI-IO, 33|. Otherwise, oscillations are relevant for the CP asymmetry, 



and it is necessary to solve the kinetic equations for the singlet neutrino propagator. In 
Section [5l we explain how such a solution may be obtained when taking the expansion of 
the Universe into account and how this solution can be used to calculate the final value 
of the lepton asymmetry of the Universe. 

Within our results, we recover many elements that have previously been derived using 
Boltzmann equations and S'-matrix elements or the Hamiltonian approach. However, 
there few features that differ from the earlier approaches: 

• For Leptogenesis from the decays of Majorana fermions, we show that the kinetic 
evolution equations fl26l) can accurately reproduce the classic result BleHiolliij for 



the asymmetry in Resonant Leptogenesis fl68|l . This requires a detailed analytical 
computation of the evolution of the two helicity states of the singlet neutrinos, and 



has not been provided in earlier work based on the Hamiltonian approach |24|-|26 



Our results consistently include the quantum statistical (Fermi-Dirac and Bose- 
Einstein) corrections that appear for all on-shell particles in the finite-density back- 
ground. Note that these also apply to the on-shell particles within the loops of the 
diagrams that describe CP-violation. 

When AM ~ Fq or smaller, it is important to take account of the lepton and 
Higgs chemical potentials, because in this situation, the lepton number violation 
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by the singlet neutrinos proceeds so slowly, that we have to take account of its time- 
dependence. This may be crucial in order to treat the washout of the asymmetry 
in a quantitatively correct manner. 

Our effective theory is formulated within one single framework, the CTP formal- 
ism. In contrast, earlier approaches typically combine elements from S-matrices, 
classical Boltzmann equations or Hamiltonian transition amplitudes. In particu- 



lar, in combination with Ref. [3J|, the methods presented here cover both, direct 
CP-violation and CP-violation from mixing, even when oscillations are important. 
Working within a single formalism may have the advantage, that it is easier to im- 
plement controlled approximations. The CTP approach combines the advantage of 
the S'-matrix approach, that perturbative expansions can be efficiently formulated 
in terms of Feynman diagrams, with the advantage of the Hamiltonian real-time 
evolution, that avoids unitarity violation from the overcounting of real intermediate 
states. 

Since the CTP approach is formulated in terms of Green functions and Feynman 
diagrams, this method for calculating the lepton asymmetry may be particularly 
suitable to be combined with corrections due to scattering processes in the thermal 



bath 38-44 . 



The CTP approach appears as a well-suited method for a consistent description of 
macroscopic statistical systems that consist of High Energy Particle Physics degrees of 
freedom. In particular, it may provide an intuitive and accurate framework in order to 
describe the emergence of the baryon asymmetry of the Universe from the quantum effect 
of CP-violation. Further significant progress on the theory of Leptogenesis appears to 
be a realistic prospect. For Resonant Leptogenesis, the present work is intended to pave 
the way to improve the quantitative predictions and the qualitative understanding. 

Additional Notes 

Another successful calculation of the asymmetry for Resonant Leptogenesis in the CTP 
approach, using partly different calculational strategies than those in the present work, 



was reported by Garny et al. [72[. During the final stage of the preparation of this 
manuscript, Ref. 73| appeared. Based on the Hamiltonian approach, it is also pointed 
out there, that the lepton chemical potential is relevant for computing the correlations 
of the singlet neutrino fiavours. 
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A Two-Point Functions on the CTP 



Two-point functions on the CTP are endowed with two path indices +/— . We often 
make use of the definitions 

G+^{x,y) = G'^{x,y), (103a) 

G+-{x,y) =G<{x,y) , (103b) 

G-+{x,y)=G>{x,y), (103c) 

G"-(x,y)=G^(x,2/), (103d) 

where G stands for any two-point function, i.e. fermionic or scalar, or Green function 
or propagator. 

The propagators of a scalar field x ^-re 

i^^{x,y) = {T[x{x)x\y)]). (104a) 

i^<{u,v) = {x\y)x{x)), (104b) 

i^>{x,y) = {x{x)x\y)), (104c) 

iA^(x,i/) = (T[x(x)xHy)]). (104d) 

In order to describe mixing of flavours, x ought to be considered as a column vector 
in flavour space. T (T) stands for (anti-)time ordering and the superscript for the 
corresponding Green functions. We sometimes refer to the Green functions with the 
superscripts <, > as Wightman functions. When x is real, the neutrality constraint 
X'^ = = X must be observed. It then follows that 

A(x,y)=A*(y,x), (105) 

where the superscript t denotes the transposition which acts here on both, CTP (±) 
and flavour indices. For both, real and complex fields the Wightman functions have the 
hermiticity property 

iA<'>(x, y) = (iA<'>(y, x)f . (106) 
Similarly, for a fermion ip, the Green functions are 

iS^{x,y) = {T[i;{x),piy)]), (107a) 

iS<{x,y) = -{i>{y)ij{x)), (107b) 

iS>{x,y) = {ij{x)^{y)), (107c) 

iS^{x,y) = {T[i:{xMy)]). (107d) 
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Notice that these define matrices in Spinor space. A Majorana spinor mus observe the 
condition 

iIj'^ = Ci}' = , (108) 

where C is the charge conjugation matrix. As a consequence, for Majorana two-point 
functions, 

S{x,y) = CS\y,x)C\ (109) 

where the transposition acts on all, CTP (±), spinor and fiavour indices. The hermiticity 
property of the Wightman function reads now 

i7°5<'>(x,2/)= (i7°5<'>(y,x))^ (110) 

which holds for both, charged and Majorana spinors. 

Additional useful two-point functions follow when taking the combinations 

G^ = G^ -G> = G<- G^ 

G^ = G^ -G< = G> - G^ 

6-^ = i(G« + G^) 

G^ = -(G^ - G^) = -(G> - G<) 
2i^ ' 2^ ' 

The Wigner transform of a two-point function is defined as 

G{k,x) = j d^re'^'G{x + r/2,x-r/2), (112) 

where we refer to r as the relative and to x as the average coordinate. In the situation of 
spatial homogeneity, there is no dependence on x (spatial translation invariance), such 
that G{k,x) = G{k,t). 

The transformation of convolution integrals in coordinate space into Wigner space 
can be expressed with the help of the diamond operator, which is defined as 

o{A}{B} = i (d^A) {dkB) - i (d^A) (d^B) , (113) 

where A = A{k, x) and B = B{k, x). 



(advanced) , 


(111a) 


(retarded) , 


(111b) 


(Hermitian) , 


(111c) 


(anti-Hermitian, spectral) . 


(Hid) 



B Equilibrium Self-Energies 



The spectral self-energy, that is central in calculations of the resonantly generated lepton 
asymmetry is computed in Ref. 33|]. Here, we generalise that expression in order to 
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account for the lepton and Higgs chemical potentials. This may be of importance when 
AM ^ Fd is not valid, such that lepton number violation is modulated by the slow 
oscillations of the right-handed neutrinos. Neglecting tree-level and thermal masses for 
the lepton and the Higgs boson, the normalised spectral self-energy can be expressed as 



NL,R 



^iVL.R(^) 



X 



167r|k| 



167r|k| 



\k^\ |k| 



rj-] ' rj-] ' I rj-\ ' i rj-\ 



|k| , 



T ' T 



T 



T 
(114a) 



(114b) 



-h 



-h 



k°\ |k| fij fi^ 



2|k|T 
2|k|T' 



W /f^ 



|fc°| |k| 



where 



h(yo+y) 



dxx'' 1- 



+ 



1 



(115) 



These integrals are given by 



Io{yo,y)= log(e^-e^O-log e^ 



(116) 



h{yo, y) = [x (log (1 + e^-^^) - log (l - e^-^^+^^j 



+ Li2 (-e"-^^) - Li2 (e 



x=Uy°+y) 

J x=\{yO-y) 



where Li2 is the dilogarithm. Of course, for ni = jj^ = 0, S'^l = ^^r- 

For the scalar model, we do not use the explicit form of the equilibrium self-energies 
for the studies presented in this paper. For completeness, we note that these can be 
calculated by defining 

pO = -i- |a;0(A;2 + ml- ml) ± •^k2(P + ml - mlf - A]^^k^ml^ (117) 
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and 

P{k,G,ma.mb,Sa,Sb,Lf) 
= G|^sign^^ (A:2) J ^{-e)k'> - pI + P° + ^ log 



/3 



log 



(118) 
e/5(fc°-P^) _ sb 



Sb 



for >0. 

The equilibrium spectral self-energy for the neutral scalar field x is then given by 



XiA I I ' 167r 



,m^,m<^,+l,+l,0 1 . 



(119) 



C Green Functions for Leptogenesis 

The propagators for the leptons are given by 

iS<{p) = -27r6{p')P^^Pn [^{po)Mp) - ^(-Po)(l - M-p))] , 
-27r5{p')P^^Pn [-^{po){l - fi{p)) + n-Po)fi{-p)] , 



i|^ 



iSjip) = Pl^^Pr - 27r6{p')P^^Pn [^{Po)fi{p) + ^{-Po)fi{-p 

iSjip) = -P^^^P^ - 27r6{p')Pi^^Pn [HPo)fi{p) + n-Po)fii-p 
p"^ — le 



(120a) 
(120b) 

(120c) 
(120d) 



Due to the unbroken SU(2)l symmetry, we suppress the indices representing the lepton 
doublet. In principle, the propagators are diagonal matrices 

:'SU{2)(^^ ^) ^ S^^S,{u, v), A,B = 1,2. (121) 



S 



Whenever these propagators appear within the integrals for the source for the asym- 
metry 5* or the singlet neutrino self-energy it is a sufficient approximation to take for 
fe{p) and fe{p) the equilibrium Fermi-Dirac distribution with zero chemical potential. 
Only in the washout term W, Eq. flW|) . one should correct this by introducing a lepton 
chemical potential, in order to describe a charged distribution in kinetic equilibrium cf. 
Eq. dHED. 

The propagators for the Higgs fields are 



iA<(p) = 27r6{p') [nPo)UiP) + ^(-Po)(l + Ui-P))] 
iA>(p) = 27t6{p') [^{po){l + Uip)) + ^{-Po)Ui-p 

iAj(p) 
iAl(p) 



+ 27r6{p') [&{po)UiP) + H-Po)U{-p 



P'^ + le 
i 



p^ — 15 



- + 27r6{p') [^{po)UiP) + ^(-Po)^(-P)] 



(122a) 
(122b) 

(122c) 
(122d) 
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As for the leptons, it is understood that 

A^V' ^) = ^AB^^u, v) A,B = 1,2. (123) 

Again, we can substitute equihbrium distributions in the colhsion integrals. When taking 
account of spectator effects, one should also assign chemical potentials to the Higgs field 
and the remaining Standard Model degrees of freedom. 

When assuming thermal equilibrium distributions and vanishing chemical potentials, 
the Kubo-Martin-Schwinger (KMS) relations are 

i^>(A;) = -e^-V^i^<(fc) , (124a) 

iA>ik)=e''/^iA<ik). (124b) 

Notice that the according relations also hold for the fermionic and scalar equilibrium 
self-energies. 



D Self-energies from the 2PI Effective Action 

The two-loop contribution 2PI effective action arising from the Yukawa interaction in 
Eq. (jH]) is given by 



= -gu,Y;Y, 5^ a6 /" rfWytr [Pr^J'(x, y)PL^^(2/, x)] A'^iy^x) , (125) 

a,b=± 

where a, 6 = ± denote the CTP indices, the trace is over Dirac indices and sum over 
neutrino flavours i is understood. The one-loop neutrino self-energy is then given by 



=g^ iY^^*P^iSfix, y)PKiAf{x, y) + Y:Y,C P^iS\''\y , x)C^P^iAl\y, x) 

(127) 

Two contributions arise from the functional derivative, because S""^ and S^"" are not in- 
dependent but related by the Majorana symmetry: Sfj{x,y) = CSj°'{y,x)C'' . Similarly, 

using gwSi = SabSi^^^ in Eq. (11251) we get for the the lepton self-energy 

=^b—^^j^^-— = Y:Y,P^iS1^MP^\A',\v,u)5ab 
=iEf{x,y)6AB. (128) 
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E Explicit Form of the Helicity Block- Diagonal De- 
composition 

In order quickly see the relation of the functions Qah with the vector, scalar, pseudoscalar 
and pseudovector densities, it is useful to recast Eq. (156|) as 



i^ 
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h=± 
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il 
-il 

-1 
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9ih 



hk' 
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-a' 
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(129) 
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